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Abstract 

In this poster I briefly review several articles on astrophysics of old isolated neutron stars, 
which were published in 1994-99 by my co-authors and myself. 



Simple Model of 


Spin And Its 


RXJ0720.4-3125 


Accretion Onto 


Evolution For 


As A Possible 


Isolated Neutron 


Isolated Neutron 


Example Of The 


Star With 


Stars: The Splndown 


Magnetic Field 


Synchrotron 
Cooling. 


Theorem. 


S>ecay In 
Neutron Stars. 


(Popov 1S34) 


(Llpunov, 
Popov 1335) 


(Konenkov, 
Popov 1337, 1338} 


2 


3 


4 




/• 


Spatial Distribution 


ABSTRACT 


Nature Of The 


Of the Accretion 




Compact X-ray 


Luminosity Of 




Source In SNR 


Isolated Neutron 




RCW 103. 


Stars And Black Holes. ^ 


- 


^ 


(Popov, 

Prokhorov 1936, 19Si 


J! 


(Popov 1338) 


5 


1 


6 


/ s, \, 


Population Synthesis 


Population Synthesis 


ROSAT X-ray 


Of Old Neutron Stars 


Of Old Neutron Stars 


sources And 


In The Galaxy: 


In the Galaxy: 


Exponential Field 


No Field Decay. 


filed Decay. 


Decay In Isolated 
Neutron Stars. 


(Popov, Colpl, 
Treves, Turolla, 

Llpunov, Prokhrov 
1333) 


(Popov, Colpl, 
Treves, Turolla, 
Llpunov, Prokhorov 
1333) 


(Popov, Prokhorov 
1333) 


7 


8 


9 



Figure 1: Plan of the poster 
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Simple Model Of Accretion Onto Isolated Neutron Star With Synchrotron Cooling 
Popov S.B., Astron. Circ. N1556 pp. 1-2, 1994 

Here I modeled accretion onto an isolated neutron star (INS) from the interstellar medium in the case of spherical 
symmetry for different values of the magnetic field strength, ambient gas density and NS's mass. I tried to verify the idea 
that if the radius of corotation, R^o, is less than the Alfven radius, Ra, the shell will form around the INS and Ra will 
decrease to Rcoi and the periodic X-ray source will appear (see Treves et al., 1993 A&A 269, 319). 

Dependence of Ra on t in our model roughly coincides with the analytic formula from Treves et al. (1993): R ^ t^^/^ 
(for some values Ra was decreasing faster). 
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Figure 2: Growth of density and decreasing of the Alfven radius 



On the figure I show the growth of the envelope density (curves on the figure arc plotted for different moments of 
time: higher density corresponds to later moments of time) and the decreasing of the Alfven radius with time. 
Periodic sources with P from several minutes to several months can appear. 



Spin and its evolution for isolated neutron stars: the Spindown theorem 
Lipunov V.M. &: Popov S.B., AZh 72, N5, pp. 711-716, 1995 ( |astro-ph/9609185| ) 

A possible scenario of spin evolution of isolated neutrons stars is considered. 

The new points of our consideration are (all points, including the Spindown theorem, are formulated for constant 
field!): 

-we give additional arguments for the relatively short time scale of the Ejector stage ( « 10^ — 10^ yrs for small velocities 
ofNSs). 

-we propose specific SPINDOWN THEOREM and give some arguments for its validity. This theorem argues, that the 
Propeller stage is always shorter than the Ejector stage (for constant magnetic field). 

-we consider evolution of spin period of a NS on the Accretor stage and predict that its period without field decay is 
> 5 ■ 10^ sec and INSs can be observed as pulsating X-ray sources. 

-we consider new idea of stochastic acceleration of very old NSs due to accretion of turbulizated ISM. A specific equilibrium 
period can be reached. 
- accreting INSs can be spin-up and spin-down with equal probability. 
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Figure S: P — y diagram 



The figure illustrates magneto-rotational evolution of an isolated neutron star on P — y diagram. Gravimagnetic 

M 



parameter: y — ^. 



RX J0720.4 3125 as a Possible Example of the Magnetic Field Decay of Neutron Stars 
Konenkov D.Yu. &: Popov S.B., PAZh 23, pp. 569-575, 1997 ( |astro-ph/9707318| ) 



Popov S.B. & Konenkov D.Yu., Radiofizika 41, pp. 28-35, 1998 (astro-ph/9812482) 



We studied possible evolution of the spin period and the magnetic field of the X-ray source RX J0720.4-3f 25 assuming 
this source to be an isolated neutron star accreting interstellar medium. Magnetic field of the source is estimated to be 
10^ — 10^ G, and it is difficult to explain observed spin period 8.38 s without invoking hypothesis of the magnetic field decay. 
We used the model of ohmic decay of the crustal magnetic field. The estimates of accretion rate (10""'^'* — lO~^^M0/yr), 
velocity of the source relative to interstellar medium (10 — 50 km/s), neutron star age (2 • 10^ — lO""^" yrs) are obtained. 

We also make new estimate of the equilibrium period for accreting INS in the ISM: 
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The period P^q corresponds to the NS rms rotation rate obtained from the solution of the corresponding Fokker-Planck 
equation. In reality, the rotational period of INS fluctuates around this value. We take into account the three-dimensional 
character of turbulence, i.e. the fact that the vortex can be oriented not only in the equatorial plane but also at any angle 
to this plane. In this case, diffusion occurs in the three-dimensional space of angular velocities. 
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Figure 4: The evolutionary tracks of the neutron star for the accretion rates M = 10^^^ Mq yr~^ (a) and M — IO^^^Mq 
yr~^ (b). The dashed lines correspond to p = Pe; the dot-dashed fines, to p = Pa- The dashed line in the figure shows 
for the second track the neutron star evolution with no acceleration in the turbulized interstellar medium. The numbers 
near the marks in tracks denote the logarithm of the neutron star age in years. The observed radio pulsars are indicated 
by dots. 



Spatial distribution of the accretion luminosity of isolated neutron stars and black holes in the GalsLxy 
Popov S.B. & Prokhorov M.E., A&A 331, pp. 535-540, 1998 ( |astro-ph/970523e| ) 



Popov S.B. & Prokhorov M.E., astro-ph/9606126, 1996 



We present here a computer model of the spatial distribution of the luminosity, produced by old isolated neutron stars 
(NS) and black holes (BH) accreting from the interstellar medium. 

We solved numerically the system of differential equations of motions in the Galactic potential. The density in our 
model is constant in time. In our model we assumed that the birthrate of NS and BH is proportional to the square of the 
local density. Stars were assumed to be born in the Galactic plane (Z=0) with circular velocities plus additional isotropic 
kick velocities. Kick velocities were taken equal for NS and BH. It is possible however that BH have lower kick velocities 
because of their higher masses. 

We used masses M^^s = IAMq for NS and Mbh = lOAf© for BH. The radii, Rub, where the energy is liberated, was 
assumed to be equal to 10 km for NS and 90 km (i.e. 3i?g, Rg ^ 2GM/c^) for BH. 

For each star we computed the exact trajectory and the accretion luminosity. The accretion luminosity was calculated 
using Bondi's formula. Calculations used a grid with a cell size 100 pc in the R-direction and 10 pc in the Z-direction. 

As expected, BH give higher luminosity than NS, as they are more massive. But if the total number of BH is 
significantly lower than the number of NS, their contribution to the luminosity can be less than the contribution of NS. 
The total accretion luminosity of the Galaxy for Nms = 10^ and Nbh = 10^ is about 10'^^ — W^^erg/s. For a characteristic 
velocity of 200 km/s the maximum of the distribution is situated approximately at 5.0 kpc for NS and at 4.8 kpc for BH. 
For NS with a characteristic velocity of 400 km/s maximum is located at 5.5 kpc, and for BH at 5.0 kpc. This result is 
also expected because of high masses of BH. 

The toroidal structure of the luminosity distribution of NS and BH is an interesting and important feature of the 
Galactic potential. As one can expect, for low characteristic kick velocities and for BH we have a higher luminosity. 

As we made very general assumptions, we argue, that such a distribution is not unique for our Galaxy, and all spiral 
galaxies can have such a distribution of the accretion luminosity, associated with accreting NS and BH. 
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Figure 5: The accretion luminosity distribution in the R~Z plane for neutron stars for a characteristic kick velocity 200 
km/s. The luminosity is in ergs per second per cubic parsec. N^s = 10^- 



Nature of the compact X-ray source in supernova remnant RCW103 

Popov S.B., Astron. Astroph. Trans. 17, pp. 35-40, 1998 ( |astro-ph/9708044 |astro-ph/9806354| ) 



Here I briefly discuss the nature of the compact X-ray source in the center of the supernova remnant RCW 103. 
Several models, based on the accretion onto a compact object such as a neutron star or a black hole (isolated or binary), 
were analyzed. 

I showed that it is more likely that the central X-ray source is an accreting neutron star than an accreting black hole. 
I also argue that models of a disrupted binary system consisting of an old accreting neutron star and a new one observed 
as a 69-ms X-ray and radio pulsar are most favored. 
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Figure 6: Possible values of the magnetic field, B, and period, p, for the accreting NS. The dotted line (1) corresponds to 
the equilibrium period, P^q, while solid line (2) corresponds to the accretor period. Pa- 



Population synthesis of old neutron stars in the Galaxy: no field decay 



Popov S.B., Col pi M. 
( |astro-ph/9910114D 



Treves A., Turolla R., Lipunov V.M., Prokhorov M.E., ApJ 530 (20 Feb), 2000 



Isolated neutron stars (NSs) are expected to be as many as 10^-10^, ~ 1% of the total stellar content of the Galaxy. 
Young NSs, active as pulsars, comprise only a tiny fraction (^ 10^"^ — 10^*) of the entire population, and about 1,000 
have been detected in radio surveys. 

The paucity of old isolated accreting neutron stars in ROSAT observations is used to derive a lower limit on the mean 
velocity of neutron stars at birth. The secular evolution of the population is simulated following the paths of a statistical 
sample of stars for different values of the initial kick velocity, drawn from an isotropic Gaussian distribution with mean 
velocity < (V) < 550 kms^^. The spin-down, induced by dipole losses and the interaction with the ambient medium, 
is tracked together with the dynamical evolution in the Galactic potential, allowing for the determination of the fraction 
of stars which are, at present, in each of the four possible stages: Ejector, Propeller, Accretor, and Georotator. Taking 
from the ROSAT All Sky Survey an upper limit of ^ 10 accreting neutron stars within ~ 140 pc from the Sun, we infer 
a lower bound for the mean kick velocity, < V >> 200 — 300 kms~^. 

Present results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to 
-1%. 
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Figure 7: Fractions of NSs in the different stages vs. the mean kick velocity for /Z30 = 0.5 (opaque circles) and fi3o — 1 
(filled circles); typical statistical uncertainty for ejectors and accretors is ~ 1-2%. 



Population synthesis of old neutron stars in the Galaxy: with field decay 



Popov S.B., Col pi M. 
( |astro-ph/9910114D 



Treves A., Turolla R., Lipunov V.M., Prokhorov M.E., ApJ 530 (Feb 20), 2000 



The time evolution of the magnetic field in isolated NSs is still a very controversial issue and no firm conclusion has 
been established as yet. A strong point is that radio pulsar observations seem to rule out fast decay with typical times 
less than « 10 Myr, but this does not exclude the possibility that B decays over much longer timescales {td ~ 10^ — 10^° 

yr). 

The same conclusion as at the previous page is reached for both a constant {B ~ 10^^ G) and a magnetic field decaying 
exponentially with a timescale ~ 10^ yr. 

We refer here only to a very simplified picture of the field decay in which B{t) = B{Q) exp {—t/td). Calculations have 
been performed for td = 1.1 x 10^ yr, td = 2.2 x 10^ yr and /X3o(0) = 1. Results are shown in the figure. 

As it is expected, the number of Propellers is significantly increased with respect to the non-decaying case, while 
Ejectors are now less abundant. Georotators are still very rare. The fraction of Accretors is approximately the same 
for the two values of td, and, at least for low mean velocities, is comparable to that of the non-decaying field while, at 
larger speeds, it seems to be somehow higher. This shows that the fraction of Accretors depends to some extent on how 
the magnetic field decays. By contrast, a fast and progressive decay of B would lead to an overabundance of Accretors 
because this situation is similar to "turning off" the magnetic field, i.e., quenching any magnetospheric effect on the 
infalling matter. 

Summarizing, we can conclude that, although both the initial distribution and the subsequent evolution of the magnetic 
field strongly influences the NS census and should be accounted for, the lower bound on the average kick derived from 
ROSAT surveys is not very sensitive to _B, at least for not too extreme values of td and /i(0), within this model. 
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Figure 8: Fractions of NSs in the different stages vs. the average kick velocity for a decaying field with an e-folding time 
td — 2.2 X 10^ yrs (opaque circles) and td — 1.1 x 10^ yrs (filled circles). 



ROSAT X-ray sources and exponential field decay in isolated neutron stars 



Popov S.B. & Prokhorov M.E., astro-ph/9908212, 1999 



Many astrophysical manifestations of neutron stars (NSs) are determined by their periods and magnetic fields. Mag- 
netic field decay in NSs is a matter of controversy. Many models of the magnetric field decay have been proposed. 

The influence of exponential magnetic field decay on the spin evolution of isolated neutron stars is studied. The 
ROSAT observations of several X-ray sources, which can be accreting old isolated neutron stars, are used to constrain 
the exponential decay parameters. Even if all modern candidates are not accreting objects, the possibility of limitations 
of magnetic field decay models based on future observations of isolated accreting NSs should be addressed. 

We show that the range of minimum value of magnetic moment, /if,, and the characteristic decay time, td, ^ 10^^'^ > 
jj^b > lO^'* G cm'^ , ^ lO'* > td > 10^ yrs are excluded assuming the standard initial magnetic momentum, /io = 10"^" G cm^. 
For these parameters, neutron stars would never reach the stage of accretion from the interstellar medium even for a low 
space velocity of the stars and a density of the ambient plasma. The range of excluded parameters increases for lower 
values of ^o ■ 

In fact the limits obtained are even stronger than they could be in nature, because we did not take into account that 
NSs can spend some significant time (in the case with field decay) at the propeller stage. 

We conclude that the existence of several old isolated accreting NSs observed by ROSAT (if it is the correct interpre- 
tation of observations), can put important bounds on the models of the magnetic field decay for isolated NSs (without 
influence of accretion, which can stimulate field decay). These models should explain the fact of observations of ~ 10 
accreting isolated NSs in the solar vicinity. Here we can not fully discuss the relations between decay parameters and 
X-ray observations of isolated NSs without detailed calculations. What we showed is that this connection should be taken 
into account and made some illustrations of it, and future investigations in that field would be desireble. 



10" 



10^ 



CO 



E 



CO 

o 

0) 
Q 



10" 



10' 



10^' 




10 10'" 10" 

Bottom magnetic momentum 



Figure 9: Characteristic time scale of the magnetic field decay, td, vs. bottom magnetic momentum, fib- In the filled region 
tE is greater than 10^°yrs. Dashed line corresponds to tn =td- In (fio/fib), where tn = 10^° years. Solid hne corresponds 
to psilJ-b) = p{t = tcr), where tcr = td ■ ln(/io/Mfc)- Both lines and filled region are plotted for /zq = 10^°Gcm~^. Dot- 
dashed line is the same as the dashed one, but for /ig = 5 • 10^^ Gcni'^. Dotted line is a border of the 'forbidden" region 
for /io = 5-1029Gcm3. 



